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Abstract: 



In a recent paper, Wisdom (2007, Icarus, in press) derived concise expressions for the rate of 
tidal dissipation in a synchronously rotating body for arbitrary orbital eccentricity and obliquity. 
He provided numerical evidence than the derived rate is always larger than in an asymptotic non- 
synchronous rotation state at any obliquity and eccentricity. Here, I present a simple mathematical 
proof of this conclusion and show that this result still holds for any spin-orbit resonance. 
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1. Introduction 

Tidal heating is a source of energy which can have a strong influence on the thermal and 
internal history of celestial bodies. For solid bodies like the terrestrial planets or rocky satellites, 
local heating is presumed to arise from the conversion of the mechanical strain energy associated 
with time-dependent tidal distortion. This may occur when a rotating satellite librates on an 
eccentric orbit or has a non-zero obliquity. 

For a homogeneous and incompressible synchronously rotating satellite, the expression of tidal 
heating as a function of the eccentricity has been calculated in detail by Peale and Cassen (1978). 
The derivation is made by calculating the power dissipated by the tide-raising force on each internal 
displaced constituent of the satellite. More recently. Wisdom (2004) generalized this expression to 
the second order in obliquity. 

For giant gaseous planets that are not expected to be trapped in spin-orbit resonances but 
rather to reach an asymptotic nonsynchronous state, Levrard et al.(2007) argued that previous 
expressions are inadequate in this situation and derived new expressions of tidal dissipation for an 
asymptotic nonsynchronous rotation valid at arbitrary eccentricity and obliquity. Unfortunately, 
they compared their results to the expressions given in Wisdom (2004) and concluded that the 
rate of tidal heating in synchronous rotation is always lower than in an asymptotic nonsynchronous 
state. Using new derivations and useful formulae for the rate of tidal heating in a synchronously 
rotating body valid for any eccentricity and obliquity. Wisdom (2007) found the opposite result. 

In this very short note, I provide a mathematical proof of this conclusion. Because solid 
(exo) planets and satellites are also expected eventually to despin to a more general state of spin- 
orbit resonance where the orbital period is some integer or half-integer times the rotation period 
(e.g. Goldreich and Peale 1966; Dobrovolskis 2007), I generalize this result to all the other spin-orbit 
resonances. 

2. Comparison of tidal dissipation between a synchronously and asymptotically 

nonsynchronous rotating body 

Wc consider the gravitational tides raised by an host planet on a satellite (the demonstration 
also holds for a planet around its central star). We use the simplest model of tidal response, 
generally called "viscous" model as described in Mignard (1980). He assumed a constant time lag 
for any frequency component of the tidal perturbation In other words, the tidally deformed surface 
of the satellite always assumes the equipotential surface it would have formed a constant time lag 
At ago, in the absence of dissipation. In this case, the ratio l/Q where Q is the satellite effective 
tidal dissipation factor is proportional to the frequency of the tides. 

Here, I calculate the rate of tidal dissipation from the variation of the mechanical energy (rota- 
tional -I- orbital) of the satellite caused by tidallly-driven perturbations in the satellite's rotational 
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and orbital parameters (e.g. Hut 1981). In that case, the total energy is 

-4--^> a, 

where lo is the satellite's rotation rate, C is its polar moment of inertia, Mp the mass of the host 
planet (primary body), Ms the satellite mass, and a is the orbit semimajor axis. 

For a satellite locked into a synchronous resonance (1:1), it is necessary to add to the right 
hand side of equation ([T]), the external gravitational potential of the deformed satellite caused by 
its permanent quadrupole moment. Averaged over an orbital period, it is classically given by0 

V = -^H{l,e){B - A)n^ cos2-i (2) 

where -ff (1, e) is the Hansen's coefficient for the synchronous resonance, e is the orbital eccentricity, 
n is the orbital mean motion, A and B are the satellite's equatorial moments of inertia and 7 is 
the resonant angle with 

dj/dt = uj — n . (3) 

The gravitational restoring torque exerted by the planet on the quadrupole moment maintains the 
spin in the resonance. 



Combining equations ([I|),([2]) and Q, the rate of tidal dissipation within the satellite then is 

■ dE ^ duj ^ , y dj GMsM„ da 

Eudai = —jr = -Cu;x — -- H{1, e) {B - A)r? -I sm 27 - x - . (4) 



Using the equation for the rotational motion of the satellite averaged over an orbital period (e.g. 
Murray & Dermott 1999) : 

= = H{1, e) {B - A) sin 27 + Tudai (5) 
where Vudai is the mean tidal torque acting to brake the spin of the satellite, equation dH) becomes 

Etidal = -UJ T tidal + -H{l,e){B - A)n^ sm 27 x — . (6) 

2 Za^ dt 

For most of satellites and planets of the solar system, the libration period is much longer than 
the orbital period but much shorter than the typical despinning timescale so that (sin(27)) = over 
a libration period. Averaging the equation ([6]) over a libration period or over secular timescales 
longer than the libration period then leads to 

i? \ IT- \ GMsMp / da\ 

Etidal I — —UJ tidal) ^-^2 ^ \'dt / ' 



^For simplicity, a zero obliquity is assumed but the same conclusion holds at any obliquity. 
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For the "viscous" tidal model, the average tidal torque is given by (e.g. Hut 1981; Neron de 
Surgy & Laskar 1997): 



tidal) 



K 

n L 



[l + x^)me)- -2xN{e) 
n 



(8) 



where e is the satellite's obliquity (the angle between the satellite equatorial and orbital planes) 
X = cos £ and 



\ I ^ I i I n (9) 

where /c2 is the potential Love number of degree 2, Rg is the satellite radius and Qn = (n At)^^ is 
the annual tidal quality factor. The eccentricity-dependent functions A^(e) and Q(e) are then 



0(e) 



1 + 3e2 + 
(1 - e2)9/2 



and 



N{e) 



1 + f e2 + f + 



(l-e2)6 

In the same way, we have (e.g. Hut 1981; Neron de Surgy & Laskar 1997): 



K 



GM,M„ 



N{e)x--Na{e) 
n 



(10) 



where 



Na{e) 



1 _L 31 „2 , 255p4 , 185 „6 , 25^8 



(1 _ e2)15/2 



Substituting equations ([8]) and ([TO]) into the equation ([7]) provides the average rate of energy 
dissipation 

Na{e) + i^5^(e) - 2xN{e) 



tidal 



2K 



(11) 



using a; ~ n and valid to any order in eccentricity and obliquity. 



We checked that this expression fully agrees with the equation (30) of Wisdom (2007) cal- 
culated for a homogeneous, incompressible and small and/or rigid enough body that the radial 
displacement Love number /12 is 5^2/3- Note that our derivation does not require such an hy- 
pothesis and all the uncertainties in the radial distribution of material and its physical proper- 
ties (e.g. density, compressibility, elasticity) are lumped into the k2 parameter. For a homoge- 
neous and incompressible body, the Love number of degree 2 is given by the well-known formula 
^2 = (3/2) / (1 + 19;u/(2 pgRg)) where p is the density and p. is the elastic shear modulus. For real 
rocky material, the effect of compressibility can not be actually neglected because p is close to A 
where A is classically the Lame's parameter which is a measure of compressibility (A — > 00 for an 
incompressible material). However, if deviations from incompressibility and homogeneous density 
requires only small corrections, a varying elastic shear modulus could have important consequences 
(e.g. Peale and Cassen 1978). For a homogeneous and incompressible fluid planet, k2 = 3/2 but is 
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about one order of magnitude smaller for realistic profiles of density. Futhermore, when the strong 
effect of compressibility is taken account, it is worthy to note that Love numbers could experience 
dramatic and unexpected variations (e.g. Hurford et al. 2002). 

Let us now consider the same satellite (or planet) that is not locked into a synchronous res- 
onance. This may occur because the body is essentially fluid or near-fluid and does not have a 
permanent quadrupole moment, or because its eccentricity is too large to allow a capture into a 
synchronous resonance (e.g. Goldreich and Peale 1966). 

Because the satellite's orbital angular momentum is generally much larger than its spin angular 
momentum, gravitational tides affect the spin's properties (rotation, obliquity) much faster than 
the orbit's properties. For the "viscous" model, tides ultimately reduce the obliquity to zero on 
the same time scale as the despinning (e.g. Hut 1981). Here, we consider that other separate 
mechanisms may maintain or excite the obliquity to a non-zero value (due for exemple to a capture 
in a Cassini state, to perturbations by a companion, ...). Setting {Tuiiai) = in the equation ([8]), 
the spin evolves to its asymptotic equilibrium rate of rotation 

N(e) 2x , , 

while the eccentricity and the obliquity are assumed to vary more slowly. Note that a non-zero 
eccentricity favors a super synchronous rotation (that is, a rotational period shorter than the or- 
bital period) while a non-zero obliquity favors a subsynchronous rotation. As a consequence, the 
asymptotic rate of rotation can be higher or lower than the synchronous rotation rate depending 
on the eccentricity and obliquity values. 

Once the spin reaches its "pseudo-equilibrium", tidal energy is then dissipated within the 
satellite only at the expense of the orbital energy. From equation d?]) , it comes 

(^i^daz) = -(GMpM,)/(2a2) x (da/dt) . 

The derived rate of tidal dissipation is then (see Levrard et al. 2007) 

iV2(e) 2a;2 



Etidal ) = 



Na{e) 



(13) 



It is now easy to compare the rate of tidal dissipation in a synchronously rotating body (Equa- 
tion (|lip ) with that in an asymptotic nonsynchronous rotation rate (Equation (|13p ) for any eccen- 
tricity and obliquity. The former is larger than the latter if 

We found that this condition is always verified because the previous equation is equivalent to 

[{l + x'^)Q{e)-2xN{e)f (15) 
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As a consequence, the dissipation rate in synchronous rotation is always equal to or larger 
than that in asymptotic rotation for any obliquity and eccentricity. Note that the equality holds 
when the asymptotic rate of rotation is synchronous, that is the satellite's eccentricity and obliquity 
verify 

N{e) 2x 

For a zero obliquity, it is possible only if the eccentricity is zero. The function N{e)/Vt{e) is 
monotonic and always increases with the eccentricity so that there is only one eccentricity and 
obliquity value that satisfies the previous condition. 



2.1. Generalization to other spin-orbit resonances 

Let us now assume that the satellite is locked into another spin-orbit resonance. We have 
LV ^ pn where p is an integer or a half- integer. Using equations ([7]), ([8]) and (jlOp . the rate of tidal 
dissipation for a resonant rotation is then 



Etidal ) — 2K 



2 l+x2. 



Na{e) +p' -^n{e) - 2pxN{e) 



(17) 



This rate of dissipation is larger than the asymptotic non-resonant rate of rotation if 



n(e)-2pxlV(e)>-^^, (18) 

Once again, this condition is always verified because this relationship is also equivalent to 

\p{\^x^)VL{e)-2xN(e)f (19) 

The equality holds when the asymptotic rate of rotation is equal to the rotation rate in the spin-orbit 
resonance that is 

w A^(e) 2x 



n r2(e) 1 -|- 

Only one eccentricity and obliquity verify this equality. 



V- (20) 



3. Conclusion 

We provide a mathematical demonstration that the rate of tidal dissipation in a synchronously 
rotating satellite is larger than that in a asymptotic nonsynchronous rotation state for any obliquity 
and eccentricity. 

We also show that this property still holds for the other spin-orbit resonances. It would be 
interesting to investigate whether these results are also valid for other tidal models. 
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